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FINAL  REPORT  FOR  PROJECT  49449-EG,  "ATOMISTIC-BASED 
MESOSCOPIC  CONSTITUTIVE  MODELS  FOR  HIGH  EXPLOSIVE 
CONSTITUENT  MATERIALS" 

FOREWORD 

This  Final  Report  pertains  to  work  performed  by  Thomas  D.  Sewell  and  postdoctoral 
co-workers  Marc  J.  Cawkwell  and  Amy  L.  Bauer  at  Los  Alamos  National  Laboratory  and 
the  Donald  L.  Thompson  Group  in  the  Department  of  Chemistry  at  the  University  of 
Missouri-Columbia  during  the  time  period  10/06-8/08  under  the  auspices  of  ARO  project 
49449-EG.  The  project  was  terminated  upon  the  departure  of  Sewell  from  the  Theoretical 
Division  at  Los  Alamos  in  August  2008  to  accept  a  faculty  position  in  the  Department  of 
Chemistry  at  the  University  of  Missouri-Columbia. 

LIST  OF  ILLUSTRATIONS 

Figure  1:  Chemical  structures  of  hexahydro- 1 ,3,5-trinitro- 1 ,3,5-triazine  (RDX),  and 
pentaerythritol  tetranitrate  (PETN). 


Figure  2:  Method  used  to  generate  shock  waves.  Three-dimensional  periodic  boundary 
conditions  were  used,  with  a  several  hundred- Angstrom  gap  introduced  at  one  end  of  the 
“long”  dimension  so  that  the  stationary  slab  on  the  left-hand  end  and  the  free  surface  on 
the  right-hand  end  of  the  periodic  simulation  cell  do  not  interact.  The  color  scheme  is 
only  intended  as  an  aid  for  the  eye. 


Figure  3:  Development  of  stacking  faults  for  RDX  shocked  at  Up  =  300  m/s  along  the 
[111]  direction.  Only  the  centers  of  mass  of  molecules  that  have  changed  numbers  of 
nearest  neighbors  relative  to  the  nominal  equilibrium  structure  are  shown.  The  original 
simulation  cell  dimensions  were  approximately  0.18  x  0.18  x  0.005  micron. 

Figure  4:  Generalized  '/-surface  enthalpy  profiles  for  the  cases  Pzz  =  0  GPa  and  Pzz  =  1 
GPa  along  the  Burgers  vector  determined  for  stacking  faults  generated  by  shock  waves 
propagating  along  [111]  in  RDX.  Left:  plots  including  estimates  of  precision.  Right: 
same,  but  with  error  bars  removed  for  ease  of  interpretation.  It  is  necessary  to  include  the 
component  of  stress  normal  to  the  stacking  fault  plane  in  order  to  render  the  stacking 
fault  metastable. 

Figure  5:  All-atom  view  of  stacking  faults  in  RDX  shocked  along  [111].  y=c/a  (lattice 
lengths). 

Figure  6:  Schematic  depiction  of  the  shock  front  absorbing  boundary  condition  approach 
used  to  the  sustain  the  initial  volume/energy  shock  state  for  times  beyond  that  required 
for  shock  wave  traversal  of  the  MD  simulation  cell.  The  essence  of  the  method  is  to 
capture  between  two  rigid  pistons  the  state  of  the  material  at  the  instant  in  time  that  the 
shock  wave  reaches  the  “free  surface”  of  a  normal  non-equilibrium  simulation  cell,  to 
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prevent  reflected  waves  from  propagating  back  into  the  cell.  (See  Refs.  2  and  3  for 
details.) 

Figure  7:  Snapshots  from  simulation  of  homogeneously  nucleated  shear  bands  in  RDX 
crystal  shocked  along  [100],  (a)  system  at  the  moment  the  shock  wave  reaches  the  right- 
hand  end  of  the  simulation  cell,  at  which  point  the  shock  front  absorbing  boundary 
condition  was  applied;  (b)  system  at  a  time  100  ps  after  application  of  the  shock  front 
absorbing  boundary  condition.  CR  and  SB  denote  regions  of  material  from  crystalline 
regions  and  shear  band  regions,  respectively,  that  were  used  for  subsequent  detailed 
analysis. 

Figure  8:  Temperature  versus  time  for  (lOO)-oriented  RDX  crystal  shocked  to  P=9.7  GPa. 
Results  are  shown  for  material  in  a  shear  band  and  for  material  in  a  crystalline  region. 
Both  intramolecular  (vibrational)  and  intermolecular  temperatures  are  shown. 

Figure  9:  Two-dimensional  map  of  temperature  behind  a  9.7  GPa  shock  in  (lOO)-oriented 
RDX  crystal.  Color  bar  units  in  Kelvin.  The  shock  front  has  passed  from  left  to  right 
through  the  material.  Defect  structures  evident  in  the  temperature  map  are 
homogeneously  nucleated  shear  bands.  Arrows  indicate  distances  of  5  nm  and  140  nm 
behind  the  shock  front.  The  original  simulation  cell  was  -0.18  x  0.18  x  0.005  micron. 

Figure  10:  Normalized  probability  distribution  functions  (PDFs)  for  temperature  in  (100)- 
oriented  RDX  crystal  shocked  to  9.7  GPa.  Results  are  shown  for  thin  slices  of  material 
perpendicular  to  the  direction  of  shock  propagation,  for  several  distances  (or, 
equivalently,  times)  behind  the  shock  front.  Top:  PDFs  for  six  distances  ranging  from  5 
nm  to  160  nm  behind  the  shock  front  (dark  blue  to  yellow,  respectively)  based  on  the 
analysis  of  a  single  simulation  snapshot  corresponding  to  t  =  36  ps.  Bottom:  as  above,  but 
with  axes  rotated  and  averaging  over  all  equivalent  time  origins  in  a  frame  centered  on 
the  moving  shock  front. 

Figure  11:  Profiles  of  intermolecular  and  intramolecular  temperature  as  functions  of 
distance  behind  the  shock  front  for  shock  waves  with  Up  =  1  km/s  propagating  along  the 
[100],  [111],  and  [210]  directions  in  RDX  crystal. 

Figure  12:  Preliminary  result  from  an  in-progress  simulation  of  shock  wave  propagation 
along  [100]  in  single-crystal  PETN.  The  piston  velocity  is  Up  =  1  km/s.  The  shock  front 
is  midway  through  the  crystal;  nucleation 
of  plasticity  at  distances  significantly 
behind  the  shock  front  is  evident. 

STATEMENT  OF  THE  PROBLEMS 
STUDIED 

The  main  goal  of  this  project  was  to 
obtain  on  the  basis  of  the  analysis  of  large- 
scale  unreactive  molecular  dynamics 
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Figure  1.  Chemical  structures  ofhexahydro-1,3,5- 
trinitro-l,3,5-triazine  (RDX),  and  pentaerythritol 
tetranitrate  (PETN). 
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(MD)  simulations  a  molecular-level  understanding  of  plasticity  and  relaxation 
mechanisms  in  polyatomic  molecular  systems,  and  to  use  this  information  to  guide  the 
formulation  of  mesoscopic  constitutive  models  for  the  constituents  of  energetic  materials 
formulations.  The  principal  focus  was  on  the  cyclic  nitramine  hexahydro-l,3,5-trinitro- 
1,3,5-triazine  (RDX);  studies  of  the  nitrate  ester  pentaerythritol  tetranitrate  (PETN)  were 
also  undertaken  near  the  end  of  the  project.  (See  Fig.  1  for  sketches  of  the  RDX  and 
PETN  molecular  structures.)  The  focus  was  on  shock  wave  loading  to  states  above,  but 
fairly  close  to,  the  Hugoniot  elastic  limit  for  oriented  crystals  of  material  without  any  pre¬ 
existing  defects;  preliminary  studies  were  performed  for  the  case  of  RDX  crystals 
containing  cylindrical  voids  in  a  quasi  2-D  simulation  cell.  (Figure  2  contains  a  depiction 
of  the  way  in  which  shock  waves  were  generated  in  these  studies.1) 

Our  approach  to  the  problem  included  three  principal  components: 

•  Identify,  for  shock  propagation  along  a  given  direction  in  a  given  crystal,  the 
dominant  mechanism(s)  of  inelastic  deformation  for  shock  loading  to  a  given 
state  on  the  Hugoniot  locus. 

•  Characterize  these  mechanisms  and  the  resulting  energy  localization  in  terms 
of  molecular-scale  phenomena  by  using  descriptors  that  capture  both  the 
“standard”  translational  degrees  of  freedom  associated  with  inelastic 
deformation  in  metals  as  well  as  additional  ones  specific  to  molecular  entities 
that  are  required  as  a  consequence  of  their  non-spherical  shapes, 
intramolecular  flexibility,  and  in  some  cases  point-group  polymorphism  or 
even  chemical  transformation. 

•  Quantify  the  practical  outcomes  of  these  inelastic  deformation  processes  on 
mesoscopic  distributions  of  thermodynamic  variables  (temperature,  stress) 
through  the  use  of  probability  distribution  functions  that  capture  the  evolution 
of  those  variables,  and  correlations  between  those  variables,  as  functions  of 
distance  or,  equivalently,  time  after  shock  passage. 

During  the  course  of  the  research  we  discovered  that  these  complicated  materials 
deform  via  pathways  that  are  in  some  ways  similar  to  those  for  metals.  For  example,  we 
observed  the  formation  of  dislocations,  stacking  faults,  and  shear  bands  depending  on 
which  material  was  studied  and  the  direction  and  intensity  of  shock  loading.  However,  we 
also  confirmed  our  expectation  that  the  observed  behaviors  include,  and  in  some  cases  are 
dominated  by,  features  for  which  there  are 
no  analogues  in  metals.  An  additional 
complication  that  had  to  be  addressed 
during  this  work  was  the  long  time  scale, 
relative  to  the  shock  wave  transit  time 
across  a  molecular  dynamics  simulation 
cell,  required  for  relaxation  of 
homogeneously  nucleated  defects  that 
developed  under  certain  loading 
conditions.  An  approach  to  dealing  with 
this  problem  was  conceived  by  Prof.  D.  L. 

Thompson  at  The  University  of  Missouri- 
Columbia  (MU)  during  the  time  frame  of 
this  project;  a  collaborative  effort  between 


Figure  2.  Method  used  to  generate  shock  waves. 
Three-dimensional  periodic  boundary  conditions 
were  used,  with  a  several  hundred- Angstrom  gap 
introduced  at  one  end  of  the  “long”  dimension  so  that 
the  stationary  slab  on  the  left-hand  end  and  the  free 
surface  on  the  right-hand  end  of  the  periodic 
simulation  cell  do  not  interact.  The  color  scheme  is 
only  intended  as  an  aid  for  the  eve. 
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MU  and  LANL  resulted  in  an  initial  method 
development  paper  followed  by  an  application 
directly  related  to  the  present  project.2,3 
Finally,  probability  distribution  functions  were 
constructed  that  capture,  on  the  basis  of 
analyses  of  simulations  containing  -500,000 
fully  flexible  molecules,  the  evolution  of 
temperature  in  RDX  as  a  function  of  distance 
behind  the  shock  front  in  that  material.  This 
kind  of  information,  and  in  particular  the 
evolution  of  the  high  temperature  tails  of  the 
distribution  corresponding  to  extreme 
fluctuations  in  the  material,  are  important  to 
the  development  of  advanced  mesoscale 
models  of  energetic  materials  initiation  under 
development  by  Mel  Baer4,5  at  Sandia 
National  Laboratory  and  Jack  Reaugh6  at  Lawrence  Livermore,  among  others. 


Figure  3.  Development  of  stacking  faults  for 
RDX  shocked  at  Up  -  300  m/s  along  the  [1 1 1] 
direction.  Only  the  centers  of  mass  of  molecules 
that  have  changed  numbers  of  nearest  neighbors 
relative  to  the  nominal  equilibrium  structure  are 
shown.  The  original  simulation  cell  dimensions 
were  approximately  0.18  x  0.18  x  0.005  micron. 


SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

Simulation  and  characterization  of  homogeneous  nucleation  of  stacking  faults  in 
(lll)-oriented  RDX.  Simulations  of  shock  loading  along  the  [111]  direction  in  RDX 
revealed  an  abrupt  onset  of  homogeneous  nucleation  of  partial  dislocation  loops  on  (001) 
with  Burgers  vector  0.16[010],  (See  Fig.  3.)  The  onset  occurs  for  shock  pressures  1.37 
GPa  <  Prh<  1-78  GPa.  The  partial  dislocation  loops  generate  stacking  faults  that  are 
obstacles  for  the  glide  of  perfect  dislocations  on  other  slip  systems,  giving  rise  to  an 
anomalous  hardening  that  does  not  involve  the  mutual  entanglement  and  storage  of 
dislocations.  Calculations  of  the  [010]  cross-section  of  the  (001)  generalized  stacking 
fault  energy  surface  as  a  function  of  applied  pressure  show  that  this  stacking  fault  is 
rendered  metastable  by  a  stress-induced  change  in  the  molecular  conformation. 
Anomalous  hardening  via  this  mechanism  is  predicted  to  result  in  high  impact  sensitivity, 
in  apparent  contrast  to  the  conclusion  one  would  obtain  upon  application  of  Dick’s7  steric 
hindrance  model  using  slip  systems  determined  from  indentation  experiments.  The  results 
of  this  paper  demonstrate:8 

•  the  importance  of  including  the  resolved  stress  on  the  stacking  fault  plane  in  order  to 
render  the  stacking  fault  metastable;  this  effect  apparently  is  unimportant  for  metals, 
due  to  the  lack  in  that  case  of  internal  degrees  of  freedom  (flexibility)  and  non- 
spherical  particle  shapes  (atoms  versus  polyatomic  molecules).  (See  Fig.  4.) 

•  the  need  for  careful  thermalization  at  each  displacement  along  the  Burgers  vector 
when  computing  the  generalized  ysurface  enthalpy  profile,  to  relax  those 
intramolecular  and  intermolecular  interactions;  in  particular,  molecules  in  the 
neighborhood  of  the  stacking  fault  undergo  changes  in  molecular  conformation  and 
orientation  in  the  crystal  (see  Fig.  5). 

•  the  need  to  exercise  care  when  applying  the  steric  hindrance  model  to  predict 
initiation  sensitivity;  in  particular,  to  assess  the  slip  systems  that  are  actually 
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Displacement  II  [010]  ( b ) 

Figure  4.  Generalized  y-surface  enthalpy  profiles 
for  the  cases  Pzz-  0  GPa  and  Pzz  -  1  GPa  along  the 
Burgers  vector  determined  for  stacking  faults 
generated  by  shock  waves  propagating  along  [111] 
in  RDX.  Left:  plots  including  estimates  of 
precision.  Right:  same,  but  with  error  bars 
removed  for  ease  of  interpretation.  It  is  necessary 
to  include  the  component  of  stress  normal  to  the 
stacking  fault  plane  in  order  to  render  the  stacking 
fault  metastable. 


Figure  5.  All-atom  view  of  stacking  faults  in 
(lll)-oriented  RDX.  y=c/a  (lattice  lengths). 


operative  under  the  regime  of  strain  rates 
of  interest  rather  than  extrapolating 


across  orders  of  magnitude  from  static  indentation  to  shock.  [PETN  appears  to  be 
something  of  a  “pathological”  case  in  this  regard,  possibly  due  to  its  small  unit  cell 
(two  molecules)  and  relatively  high  space  group  symmetry.] 
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Extension  of  the  Shock  Front  Absorbing  Boundary  Condition  and  study  of 
homogeneously  nucleated  shear  bands  in  shocked  (lOO)-oriented  RDX  single  crystal. 

The  response  of  the  energetic  molecular  crystal  RDX  to  the  propagation  of  planar  shock 
waves  normal  to  (100)  was  also  studied.3  In  this  case,  we  observed  the  nucleation  and 
spread  of  shear  bands  during  shocks  with  a  particle  velocity  of  1.0  km/s  and 
corresponding  Rankine-Hugoniot  shock  pressure  of  9.7  GPa.  These  defects  propagate 
into  the  compressed  material  at  45°  to  [100]  in  the  (010)  zone,  and  lead  to  localized 
temperature  rise  due  to  viscous  heating  in  the  shear  bands.  In  contrast  to  the  case 
discussed  above,  however,  the  shear  bands  evolve 
slowly  compared  to  the  time  scales  required  for 
the  shock  wave  to  traverse  the  simulation  cell. 

Thus,  a  recently  developed2  shock  front  absorbing 
boundary  condition  was  applied  (Fig.  6),  with 
some  adaptation,  to  the  simulation  cells  studied 
here  to  sustain  the  shock-compressed  state.3  This 
enabled  simulations  of  sufficient  duration  to  study 
the  temporal  and  structural  evolution  of  the  shock- 
induced  shear  bands  toward  a  steady-fluctuating 
state.  Analyses  of  orientational  and 
conformational  disorder  were  performed  and 
contrasted  for  regions  within  and  external  to  the 
shear  bands  (Fig.  7).  Additionally,  an  estimate  of 
the  “shock  underheating”  due  to  the  classical 
approximation  was  performed.  The  major 


Figure  6.  Schematic  depiction  of  the  shock 
front  absorbing  boundary  condition 
approach  used  to  the  sustain  the  initial 
volume/energy  shock  state  for  times  beyond 
that  required  for  shock  wave  traversal  of  the 
MD  simulation  cell.  The  essence  of  the 
method  is  to  capture  between  two  rigid 
pistons  the  state  of  the  material  at  the  instant 
in  time  that  the  shock  wave  reaches  the  “free 
surface”  of  a  normal  non-equilibrium 
simulation  cell,  to  prevent  reflected  waves 
from  propagating  back  into  the  cell.  (See 
Refs.  2  and  3  for  details.) 
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conclusions  from  this  work  are: 

•  Due  to  intense  viscous  heating  within  them, 
the  shear  bands  can  be  considered  to  be 
homogeneously  nucleated  hot  spots  (Fig.  8). 

•  Two  of  the  three  structural  determinants  that 
characterize  liquid  behavior  —  molecular 
reorientation/conformational  relaxation  and 
pair  distribution  function  g(r)  —  were  studied 
and  confirmed  that  material  in  the  shear  band 
is  indeed  liquidlike;  the  third,  self-diffusion, 
was  not  considered  due  to  the  complicated 
counterflow  boundary  conditions  associated 
with  the  shear  bands. 

•  Both  molecular  orientations  and 
intramolecular  properties  (distributions  of 
axial  and  equatorial  NO2  groups  and 
boat/chair/half-boat  six-member  ring 
conformations)  are  sensitive  indicators  of  shear  banding. 

•  For  the  9.7  GPa  shock  strength  considered  here,  the  time  constant  for  rotational 
reorientation  in  the  shear  band  is  —15  ps. 

•  Using  the  shock  front  absorbing  boundary 
condition,  it  was  practical  to  simulate  a  system 
for  times  sufficiently  long  to  observe 
nucleation,  growth,  and  eventual  quenching 
and  cooling  ( via  thermal  conduction)  of  the 
shear  bands. 

•  A  careful  analysis  based  on  equation  of  state 
considerations  was  performed  to  estimate  the 
magnitude  of  underheating  compared  to 
experiment  for  the  simulated  shock  in  regions 
far-removed  from  shear  bands.  The 
underheating  arises  due  to  differences  in 
behavior  of  the  classical  mechanical  and 
quantum  mechanical  specific  heat  of  the 
material,  the  former  corresponding  to  the 
molecular  dynamics  simulations  and  the  latter 
to  the  real  material. 

Probability  distribution  functions  of  temperature  in  RDX  crystal  shocked  along 
[100].  The  average  temperatures  attained  by  purely  hydrodynamic  shock  heating  are 
insufficient  to  initiate  significant  chemical  reactions  in  plastic-bonded  explosives,  even 
for  bulk  shock  pressures  for  which  detonation  initiation  is  known  to  occur.  Thus,  energy 
localization  processes,  whether  homogeneously  or  heterogeneously  nucleated,  are 
required  to  achieve  sufficiently  high  local  energy  density  (“hot  spots”)  for  rapid 
chemistry  to  initiate  (Fig.  9).  Moreover,  given  the  exponential  dependence  of  chemical 
reaction  rates  with  temperature,  it  is  the  extreme  fluctuations  in  local  temperature, 


Figure  8.  Temperature  versus  time  for 
(lOO)-oriented  RDX  crystal  shocked  to 
p=9.7  GPa.  Results  are  shown  for  material 
in  a  shear  band  and  for  material  in  a 
crystalline  region.  Both  intramolecular 
(vibrational)  and  intermolecular 
temperatures  are  shown. 


Figure  7.  Snapshots  from  simulation  of 
homogeneously  nucleated  shear  bands  in 
RDX  crystal  shocked  along  [100].  (a) 
system  at  the  moment  the  shock  wave 
reaches  the  right-hand  end  of  the  simulation 
cell,  at  which  point  the  shock  front 
absorbing  boundary  condition  was  applied; 
(b)  system  at  a  time  100  ps  after  application 
of  the  shock  front  absorbing  boundary 
condition.  CR  and  SB  denote  regions  of 
material  from  crystalline  regions  and  shear 
band  regions,  respectively,  that  were  used 
for  subsequent  detailed  analysis. 
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corresponding  to  the  high-temperature 
tails  of  the  associated  distribution 
functions,  which  will  dominate  the 
earliest  stages  of  initiation  chemistry. 
(Analogous  comments  apply  to  the 
importance  of  local  stress  in  the 
shocked  material.) 

In  order  to  characterize  these 
temperature  fluctuations,  real  space 
probability  distribution  functions  of 
temperature  in  (lOO)-oriented  RDX 
crystal  shocked  to  P  =  9.7  GPa 
(sufficient  to  nucleate  shear  banding  in 
the  material;  see  preceding  section) 
were  determined  as  a  function  of 
distance  or,  equivalently,  time  after 
shock  passage  through  the  material 
(Fig.  10,  top  panel).  Thus,  each  of  the 
six  curves  in  top  panel  of  Fig.  10 
corresponds  to  the  distribution  of 
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Distance  behind  shock  front  (nm) 


Figure  9.  Two-dimensional  map  of  temperature  behind 
a  9.7  GPa  shock  in  (lOO)-oriented  RDX  crystal.  Color 
bar  units  in  Kelvin.  The  shock  front  has  passed  from  left 
to  right  through  the  material.  Defect  structures  evident  in 
the  temperature  map  are  homogeneously  nucleated  shear 
bands.  Arrows  indicate  distances  of  5  nm  and  140  nm 
behind  the  shock  front.  The  original  simulation  cell  was 
-0.18  x  0.18  x  0.005  micron. 


temperatures  observed  in  a  slice  of  material  perpendicular  to  the  direction  of  shock 
propagation,  at  distances  ranging  from  5  nm  to  160  nm  behind  the  shock  front.  Even 
though  each  of  these  distribution  functions  is  based  on  the  kinetic  energies  of  several 
thousand  molecules,  the  statistics  are  rather  noisy.  However,  by  choosing  a  reference 
frame  that  is  stationary  with  respect  to  the  moving  shock  front,  it  is  possible  to  collect 
simulation  data  over  multiple  time  origins  and  therefore  improve  significantly  the 
accuracy  and  precision  of  the  distributions  (Fig.  10,  bottom  panel).  These  improved 
statistics  bring  into  relief  the  rapid  relaxation,  within  tens  of  nm  behind  the  shock  front  in 
the  present  case,  of  an  initial  tail  in  the  temperature  distribution  distances  that  extends  out 
to  -1000  K,  even  though  the  average  temperature  behind  the  shock  wave  is  only  -400  K. 

The  results  shown  here  are  limited  to  single  variable  (temperature)  probability 
distribution  functions  (PDFs)  for  slices  of  material  perpendicular  to  the  direction  of  shock 
propagation  and  at  specific  distances  behind  the  shock  front.  While  this  kind  of 
information  is  interesting,  access  to  conditional  PDFs  will  yield  considerably  more 
information  of  use  to  mesoscale  model  development.  For  example,  determining  the 
conditional  PDF  for  observing  a  temperature  T at  some  point  y  given  a  temperature  T’  at 
some  point  y  ’,  both  at  some  distance  x  behind  the  shock  front,  will  begin  to  provide 
information  concerning  spatial  correlations  among  hot  regions  in  the  material;  extension 
to  time  evolution  of  these  spatial  correlations  (by  interpolating  between  slices 
corresponding  to  different  distances  behind  the  shock)  should  be  straightforward,  and  will 
lead  to  tangible  information  concerning  three-dimensional  hot  spot  sizes  and  shapes. 
Probability  distribution  functions  analogous  to  those  described  above,  but  for  which 
stress  is  the  variable  of  interest,  as  well  as  multivariate  PDFs  involving  both  temperature 
and  stress  (and  possibly  local  strain  or  strain  deformation  gradient),  are  also  of  significant 
interest.  Finally,  an  intriguing  possibility  is  to  study  these  PDFs  in  reciprocal  space,  since 
important  structures  in  the  data  may  not  be  apparent  in  the  real-space  picture  yet  may 
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Figure  10.  Normalized  probability  distributions 
(PDFs)  for  temperature  in  (lOO)-oriented  crystal 
shocked  to  9.7  GPa.  Results  are  shown  for  thin  slices 
of  material  perpendicular  to  the  direction  of  shock 
propagation,  for  several  distances  (or,  equivalently, 
times)  behind  the  shock  front.  Top:  PDFs  for  six 
distances  ranging  from  5  nm  to  160  nm  behind  the 
shock  front  (dark  blue  to  yellow,  respectively)  based 
on  the  analysis  of  a  single  simulation  snapshot 
corresponding  to  t=36  ps.  Bottom:  as  above,  but  with 
axes  rotated  and  averaging  over  all  equivalent  time 
origins  in  a  frame  centered  on  the  moving  shock  front. 


stand  out  clearly  in  reciprocal  space. 
These  latter  possibilities  comprise  a 
research  area  of  intense  ongoing  interest. 
Additional  Studies.  Additional  studies 
that  were  begun  during  this  contract 
include  simulations  of  shock  waves  with 
Up  =  1  km/s  propagating  along  three 
different  directions  in  single  crystal  RDX 
(Fig.  11)  and  a  simulation  of  PETN 
crystal  shocked  along  [100]  at  Up  =  1 
km/s  (Fig.  12).  In  the  case  of  the  RDX 
simulations,  plots  of  intermolecular  and 
intramolecular  temperature  profiles  as 
functions  of  distance  (or,  equivalently, 
time)  behind  the  shock  front,  determined 
for  different  directions  of  shock  wave 
propagation  ([100],  [111],  and  [210]) 
reveal  significant  differences  in  the  extent 
of  initial  overheating  of  lattice  modes 
(intermolecular  temperature)  and 
distances  required  for  the  approach  to 
thermal  equilibrium  among  the 
intermolecular  and  intramolecular 
measures  of  temperature.  Understanding 
the  origins  of  these  differences  in  terms  of 
fundamental  deformation  processes  and 
dynamical  energy  transfer  mechanisms  is 
of  keen  intellectual  interest,  and  will  lead 
to  useful  insight  concerning  intrinsic 
shock  initiation  sensitivity. 

The  result  shown  for  (lOO)-oriented 
PETN  indicates  the  onset  of  plasticity. 
The  details  of  the  deformation  have  yet  to 
be  explored,  however,  in  spite  of  the 
assessing  the  steric  hindrance  model  as  an 


current  interest  in  understanding  and 
explanation  for  anisotropic  initiation  sensitivity  based  on  the  presence  or  absence  of 
plastic  slip  systems  for  different  directions  of  shock  loading.  While  neither  of  these 
investigations  reached  a  point  that  warrants  extensive  discussion,  both  have  revealed 
intriguing  behaviors  that  should  form  the  basis  of  future  research. 
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Figure  11:  Profiles  of  intermolecular  and  intramolecular  temperature  as  functions  of  distance  behind  the 
shock  front  for  shock  waves  with  Up  =  1  km/s  propagating  along  the  [100],  [111],  and  [210]  directions  in 
RDX  crystal. 


Figure  12.  Preliminary  result  from  an  in-progress  simulation  of  shock  wave  propagation  along  [100]  in 
single-crystal  PETN.  The  piston  velocity  is  Up  -  1  km/s.  The  shock  front  is  midway  through  the  crystal; 
nucleation  of  plasticity  at  distances  significantly  behind  the  shock  front  is  evident. 
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